Fighting dephasing noise with robust optimal control 
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We address the experimentally relevant problem of robust mitigation of dephasing noise acting 
on a qubit. We first present an extension of a method for representing 1/lj" noise developed by 
Kuopanportti et al. [I] to the efficient representation of arbitrary Markovian noise. We then add 
qubit control pulses to enable the design of numerically optimized, two-dimensional control functions 
with bounded amplitude, that are capable of decoupling the qubit from the dephasing effects of a 
broad variety of Markovian noise spectral densities during arbitrary one qubit quantum operations. 
We illustrate the method with development of numerically optimized control pulse sequences that 
minimize decoherence due to a combination of 1/cj and constant offset noise sources. Comparison 
with the performance of standard dynamical decoupling protocols shows that the numerically opti- 
mized pulse sequences are considerably more robust with respect to the noise offset, rendering them 
attractive for application to situations where homogeneous dephasing noise sources are accompanied 
by some extent of heterogeneous dephasing. Application to the mitigation of dephasing noise on 
spin qubits in silicon indicates that high fidelity single qubit gates are possible with current pulse 
generation technology. 

PACS numbers: 



INTRODUCTION 



Coherent control of interactions between quantum bits 
(qubits) and their environment is an essential component 
of the search for realization of reliable quantum compu- 
tation. Exercising control of coherence by application of 
carefully designed pulse sequences is a standard tool in 
NMR and ESR, where particularly powerful sets of tools 
have been developed for protection against qubit dephas- 
ing. Such coherent control provides a complementary 
approach to the protection of quantum information by 
encoding, whether with active quantum error correction 
codes or passive encoding into decoherence free ('noise- 
less') subspaces and subsystems. In the quantum infor- 
mation processing community, the application of coher- 
ent control ideas to preservation of qubit coherence has 
prompted an explosion of work in the field of dynamical 
decoupling [21 [3] which builds on the classic spin echo 
techniques from NMR. An alternative approach is the 
design of numerically optimized control pulse sequences, 
which have the advantage of flexibility and ready ap- 
plicability to both quantum memory and protection of 
arbitrary quantum gates against dephasing and leakage 
errors [THHIS]. 

In this work we consider the design of optimal pulse 
sequences for control of decoherence during single qubit 
operations when the qubit is coupled to source of Marko- 
vian noise that possesses an arbitrary noise spectrum. 
The Hamiltonian governing the evolution of the qubit is 



taken to be 

H{t) = - {a^(t)aa, + ay{t)(Ty + r]{t)a,) , (1) 

where we are working in a rotating frame so that the 
qubit energy level splitting is, on average, zero. Here, 
a^it) and ay{t) are bounded-amplitude control fields, 
while ri{t) is a classical stochastic process. In the ab- 
sence of control, this Hamiltonian generates dephasing 
of the qubit, at a rate that depends on both the am- 
plitude distribution and the temporal correlation func- 
tion of the stochastic process ri{t). Previous work in our 
group has modeled this stochastic process as a multi- 
state Markovian fluctuator, the dynamics of which may 
be tuned to mimic a particular noise source [T] The use 
of a multi-state Markovian fluctuator allows the evolu- 
tion under the Hamiltonian, Eq. (fTl) to be efficiently 
solved through use of a deterministic master equation 
E]. In the present work, we extend this approach from 
the l/w" spectral noise sources for which analytic repre- 
sentations of the fluctuator dynamics could be found [T] , 
to representation of a broad variety of Markovian spec- 
tral noise distributions by making use of convex opti- 
mization techniques. We then use gradient ascent meth- 
ods as in Ref. [T] to derive control sequences for single 
qubit operations. Following established literature con- 
vention, these pulse sequences will be generically refered 
to as "GRAPE sequences" (i.e., gradient-ascent pulse en- 
gineering sequences) 0. We focus here on two specific 
single qubit operations: i) the identity, which is equiv- 
alent to extension of qubit coherence, i.e., to quantum 



memory; ii) the Hadamard gate. We analyze the ro- 
bustness of these pulse sequences for protection against 
the combined effects of 1/w noise and zero frequency 
noise (the latter is equivalent to a resonance frequency 
error) and compare with the corresponding performance 
of standard dynamical decoupling sequences, in partic- 
ular with Carr-Purcell sequences. We find that the nu- 
merically optimized control sequences improve on the dy- 
namical decoupling sequences over a broad range of zero 
frequency noise offsets, resulting in considerably greater 
robustness in addition to improved decoherence mitiga- 
tion. Finally, we demonstrate the applicability of the 
method for current day experiments by making explicit 
application to the protection of coherence for dopant spin 
qubits in silicon using realistic estimates of spectral noise 
and control pulse capabilities. The results indicate that 
gates with errors less than 10~^ can be designed and im- 
plemented with current technology. This is well below 
current estimates of the fault-tolerant threshold [8] . 



SIMULATED NOISE MODEL 



Qubit evolution under classical noise ry(i) is simulated 
using an extension of a method first developed in [1], 
which we review and then expand upon here. We con- 
sider r](t) in Eq. (II]) to represent a multi-state Markovian 
fluctuator having N noise states. The noise state k has 
amplitude rjk and occupation probability pfc(i) at time, t. 
These noise amplitudes and occupation probabilities will 
be represented as vectors, ff and p. Transitions between 
noise states are governed by a rate matrix, F, 



^p(t) = Fp(t). 

To conserve probability, the transition rate matrix must 
satisfy ^j, F^^ = 0. This constraint implies that the vec- 
tor ps = {l/N, 1/N, . . . , l/N) is a stationary probability 
vector, and is an eigenvector of F with zero eigenvalue. 
We shall limit our study of the rate matrices to those 
satisfying the additional requirement, F = F""". This con- 
dition makes the forward and backward transition rates 
between any two noise states to be equal, enforcing time- 
reversal invariance on the fiuctuator dynamics. 

The noise source may be further characterized by its 
power-spectral density, S{lli). 
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culated through its temporal correlation function 

c{t) = {v{t)vm 
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where P{r]i{0)) = Pi{0) denotes the probability of the 
fiuctuator being in state i at time t — and P(r]j, t\rii, 0) 
the conditional probability of it being in state j at time 
t, given state i aX t = 0. Here we have chosen as the 
initial noise probabilities the stationary vector P{t) = 
Ps = {1/N,1/N, . . . ,1/N). Because F is a symmetric 
matrix, it can be diagonalized by an orthogonal matrix, 
F = ytAV, so that 

C{t) = ^7fV^e^\'\Vff^Pe^\% (3) 

where we have defined the transformed noise amplitude 
vector b = VfJ/y/N and A = diagjAi, A2, . . . , Ajv} is the 
diagonal matrix of eigenvalues of F. For convenience, 
the eigenvalues are ordered i > j ^ Xi < Xj. The cor- 
responding power-spectral density, Eq. (pi), is a sum of 
zero-mean Lorentzian distributions: 
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(4) 



Ref. [I] derived an analytic form of F and V that gen- 
erates noise with a 1/w", < a < 2, power spectrum. 
As noted there, numerical optimization may result in a 
more accurate representation. 

As with all Markovian processes, the form of Eq. ^ 
is, in accordance with Doob's theorem [5], a sum of 
Lorentzians. This form constrains the possible target 
spectra to those which are monotonically decreasing and 
which never decay faster than l/oj^. We have found that 
by proper choice of A and b, this spectrum may indeed 
be brought arbitrarily close to a given target spectrum, 
Stitj), (chosen with the above constrains in mind) over 
a finite specified range of frequencies, uj G [wmin; Wmax]- 
The choice of A and b is made by a numerical optimiza- 
tion that minimizes the deviation of Eq. (l4| from the 
target spectrum. In particular, we carry out the follow- 
ing optimization: 



minimize 
\,b 



W{oj) (^(w; A, fe) - S't(w)) dw 



The power spectrum of the multistate fluctuator is cal- 



subjectto bi > 0, Xi < 0. 

Since i) w can span many orders of magnitude, and ii) 
analytic representations of power spectra often diverge at 
a; = 0, we have incorporated here a weighting function, 
W{uj), into the usual L^ distance measure. In particu- 



lar, we have set W{uj) = 1/ui for w G [wminiWmax] and 
W{uj) = otherwise. This weight function is uniformly 
distributed in logw, preventing higher frequencies from 
dominating the integral. Restricting ourselves to the 
range A^ < is physically realistic, since positive eigen- 
values would not conserve probability and would cause 
the correlation function Eq. ^ to diverge at long times. 



The results of this optimization are the two vectors 
6opt and Aopt. Recall that the constraints on F imply 
the existence of a stationary probability vector, ps with 
eigenvalue Ai = 0. This implies that one component of 
6opt, say 6i, can be taken to be a free parameter and may 
be chosen to make the arithmetic mean oi ff = ^/NV'b 
equal to zero, guaranteeing the existence of the stationary 
solution Ps = {l/N, 1/N, ..., l/N). We note that it is 
convenient to further make a restriction to 5^ > during 
the numerical optimization, because the power spectrum 
depends only on b^. However, following the optimization, 
we may subsequently adjust the signs of all components 
6i/i so that \bi\ is as small as possible, consistent with 
the existence of the stationary probability vector. 

It now remains to construct a valid transition rate ma- 
trix r with eigenvalues given by Aopt • This is again done 
by a numerical optimization, namely 



minimize 
r 



(eigs(r) 



A, 



opt 



subject to r = r"^, T,^j > 

3 

with eigs(r) the vector of eigenvalues of T. 

Both of the above optimizations can performed very 
quickly using standard convex optimization tools [TOllllj . 
Furthermore, this numerical approach is general and may 
be readily applied to generate other spectral distribu- 
tions (subject to the constraints mentioned above). In 
the remainder of this paper we shall employ a multi- 
state Markovian fluctuator representation of l/w noise 
together with a zero frequency component rjos that de- 
scribes possible sources of heterogeneous dephasing. This 
combination is experimentally relevant to a broad range 
of physical qubits. 
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with 77os the constant noise offset responsible for zero fre- 
quency noise, e the scaling of the noise amplitude (repre- 
sented as a fraction of the maximum control amplitude) 
and Tm a constant that tunes the range of frequencies 
over which the fluctuator best approximates St- The 
resulting power spectrum for parameter set e = 10^'^, 
Tm = 1/30 and 770s = is shown for a range of finite 
frequency oj in Fig. [ll For comparison we also show the 
corresponding approximation to the 1/w power spectrum 
derived from a 32-state Markovian fluctuator with the 
analytic form of Ref. [T]. It is evident that the numer- 
ically optimized 4-state Markovian fluctuator provides 
a signiflcantly improved fit relative to the analytic ap- 
proximation, as well as a significantly greater range of 
representation. Such enhanced accuracy, together with 
the considerable increase in efficiency and greater flexibil- 
ity, illustrated here by the addition of the zero frequency 
noise component in the fit (see Eq. ([5])), render this nu- 
merical optimization approach to generation of arbitrary 
spectral noise densities extremely attractive. 
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Figure 1; (Color online) Numerically optimized noise power 
spectral density (red, dashed line) with optimization con- 
structed to match the target noise spectrum St oc l/w -I- 
rias5{iu) (blue, solid line) over two decades of frequency uj. 
Also shown is the fit obtained with the analytic representa- 
tion of Ref. [1] using 32 noise states (yellow, dotted line). 



To illustrate the efficiency and flexibility of this gen- 
eration of arbitrary noise spectral densities by numer- 
ical optimization of a multi-state Markovian fluctua- 
tor, we applied the approach to calculation of a 4-state 
Markovian fluctuator representation of a target spectrum 
St (X l/w + ■q1^6{uj). The result is 



rf=r]os + e(-0.875, 1.36, -1.36, 0.875) 



(5) 



Qubit Evolution with Noise 

The evolution of this one qubit system under classi- 
cal dephasing noise r]{t) is exactly solvable through the 
use of conditional density matrices, Pk{t), as described 
in [5] . We outline here a slightly modifled version of this 
approach. Defining pk (t) as the density operator of the 



system conditioned on the environment being in the state 
k, the total density operator of the system is given by the 
sum of the conditional density operators weighted by the 
probability of occupation of the associated noise state: 



N 



p{t)^Y.pkii)pk{t)- 



even (and be called a "pulse"), and will take zero ampli- 
tude for i odd, corresponding to a quiescent time between 
pulses. Thus, TVp is understood to mean the number of 
pulses in the control pulse sequence. Each control has 
time duration 5i and the total time for a pulse sequence 
is equal to r = ^^ 5i. With 27Vp control function values, 
the corresponding Bloch vector dynamics are given by 



We choose as the initial probability vector p(0) the sta- 
tionary probability, pg. For a single qubit, we can pa- 
rameterize the conditional density matries Pk{t) by their 
Bloch vectors, Ck{i)i 



pu{t) = -[T + Qkit) ■ a 



where T is the identity operator and a is the vector of 
Pauli spin-1/2 matrices. The resulting stochastic Liou- 
ville equation for the conditional density matrices can 
be transformed to yield the dynamics of the conditional 
Bloch vectors, which are given by 



27V„ 



dt 



Ckit) = Mfc(i)a(i) + Y,^Mt). 



(6) 



i(r)- rJlexplA^O i(0). 



i=l 



Calculating the evolution of a given initial state is then a 
matter of matrix multiplication. Because the probability 
vector p(t) = ps, all noise states are equally probable and 
the relation between Z{t) and the Bloch vector, C(i), is 
given by 



m 
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Z{t) ^ -lNZ[t) 



The inverse relation is simply ^(0) — I]^^C(0) and the 
final Bloch vector is then C(t) = f C(0), where 



Here M^ G so (3) is the generator of Bloch vector rota- 
tions, 
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The second term in Eq. ^ describes the effect of the 
noise switching on the conditional Bloch vectors. We 
have thus arrived at a set of N coupled matrix differ- 
ential equations for the evolution of the N conditional 
Bloch vectors. These can be solved by treating the set of 
Bloch vectors as a single vector composed by stacking the 
conditional Bloch vectors to get a single 37V-dimensional 
vector, Z{t) = ®fc^iCfc(0- The equation of motion for 
Z{t) can be straightforwardly derived from Eq. ([6]) and 
is given by 






N 



0MfcW 



+ V®X\ Z{t) = C{t)Zit). 



\k=l 



This is solved formally in the usual way, namely as 
Z{t) = T exp (j^ C{t')dA Z(0), where the symbol, T, 

is the usual Dyson time-ordering operator. This time- 
ordered integral becomes a time-ordered product if we 
restrict the control functions, ax{t) and ay{t), to those 
that are piecewise-constant in time. The Lindblad op- 
erators, C{t), are also then piecewise constant, taking 
values Ci for times 5i. For future convenience, we di- 
vide the control functions into 2Np subintervals, where 
subinterval i will in general take nonzero amplitude for i 



£ = —Tn ■ T]J exp (A 5ti) • XJ^. 



(7) 



Numerically Optimized Control 

The control functions must now be chosen to realize a 
desired target operation on the Bloch vector, Q — > G^. 
We choose the operator fidelity to measure the efficacy 
of these control functions. We define the fidelity function 
as |I?j . 



$GK(t),a,W] = 2(l + C(r)-GC(0)) 



1+£C(0)-GC(0) 



Note that the fidelity is a functional of the control fields, 
o-x/yit)- From the perspective of quantum information, 
no state is any more important that any other, so we 
would ideally like our pulse sequences to maximize the 
worst-case fidelity over all possible initial states. How- 
ever, the minimization over initial states to find the 
worst-case fidelity is too expensive a computation to yield 
a useful cost function. Therefore we use instead as our 
cost function the average case fidelity and compare this 
with the worst case fidelity obtained from the optimized 
pulse sequence at the end of the computation in order to 



ascertain the range of errors. Thus we 



subject to ttxit)^ + ayit)"^ < 1, 



maximize 



where the notation, (•)?, imphes an average taken over 
the surface of the Bloch sphere. The above average can 
be evaluated as 



1 1 



a>GK(t),a,(<)])^^(0) = 9 + 9 E (^.,G.fcC,(0)a(0)>f(o) 
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We have used the resuh that the average over the surface 
of a unit sphere is given by {CiCj)r ~ ^ijl'^- Finally, as 
discussed above, we demand that the optimized pulse 
sequences be insensitive to zero-frequency noise. This is 
achieved by including a constant offset value, 770s, to the 
noise vector, as in Eq. [SJ We systematically analyze the 
effect of this additional zero frequency noise by choosing 
the offset to take values within the range \t]os\ < lOe. The 
optimization problem then becomes 



maximize 

ax(t),ay(t) 



mm 

JJos 



Tr(£G"^) 



(8) 



subject to ax{tY + ayiif < 1- 



It is in general possible to find analytic gradients of 
^G[0'x(t),ay(t)] in terms of the pulse sequence parame- 
ters by straightforward methods of |13) when rjos is fixed. 
However, because the objective function has the form of 
a minimum value over some range of 770s, the objective 
function is not in general differentiable everywhere. We 
therefore use the solver [TUl E] which employs finite dif- 
ference approximations to the gradient (which may be 
undefined in certain regions). A finite difference mini- 
mization approach requires many more function evalu- 
ations than an explicit gradient calculation, greatly in- 
creasing optimization time. 

The optimization is performed by undertaking a sam- 
pling over the allowed parameter space. We begin by 
randomly selecting an initial point in the space, and ap- 
plying numerical optimization techniques to find a locally 
optimal value of the objective function. We repeat this 
process many times, each time obtaining a value for $g. 
After some fixed number of initial conditions are sampled 
(typically thousands), the pulse sequence obtaining the 
greatest value of the operator fidelity, $g, is selected as 
the optimal sequence. 

It is important to note that the dimension of the pa- 
rameter space for optimization is 4iVp. A pulse sequence 
contains A^p pulses, each of which can be characterized 
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Figure 2: Worst-case error, (1 — "3?h), over the Bloch sphere 
for a Hadamard gate generated by GRAPE as a function of 
the number of pulses in the pulse sequence up to Np = 15 for 
fixed computational efi^ort. See text for disucssion. 



by four parameters: the amplitude, duration, phase an- 
gle of the control fields, and quiescent period before the 
next pulse. To thoroughly sample the enlarging param- 
eter space, the number of initial points sampled for op- 
timization should grow as n^ ^ ^ where n^ is the number 
of initial samples taken for iVp — 1. 

An interesting point is raised by Fig. [2] which shows 
that the worst-case error of a numerically optimized 
Hadamard gate as a function of A^p. If the computa- 
tional effort is allowed to grow exponentially with N-p 
(and therefore linearly with the size of the parameter 
space), we can expect the gate error to decrease mono- 
tonically with A^p. With exponential computational re- 
sources, one is able to sample all of the parameter space 
defining the pulse sequences. And because the set of 



N, 



N, 



n pulse sequences is a strict subset of the set of 
p — n -I- 1 pulse sequences, such a search should yield 
sequences that, at the very worst, do not decrease in effi- 
cacy. However, when the number of initial sample points 
does not grow exponentially with A^p, the performance of 
GRAPE will suffer because it becomes exceedingly un- 
likely to find good optima of the objective function as 
A^p becomes large. This exponential scaling in the num- 
ber of sample points required places an upper limit on 
the largest values of A^p for which the GRAPE approach 
will be useful. In the examples studied in this work, we 
find that for gates other than the identity, this scaling 
restricts values of A'p to a maximum of 6-10. Thus in 
Fig. [2] we see that the error begins to rise after A^p — 10. 
For non-trivial gates this is not a disadvantage since it 
is in any case advantageous to implement the gate in as 
short a time as possible to avoid decoherence. The ex- 
ception to this is the identity gate, which one wishes to 
realize over as long a time as possible when implement- 
ing a quantum memory. In this case we find that good 
identity gates may be found for A'p ^ 30. 



RESULTS 

To demonstrate the power and flexibility of the im- 
proved multi-state Markovian fluctuator approach to co- 
herent control of qubit dephasing, we apply it to study- 
ing magnetic interface noise experienced by a phospho- 
rus dopant atom implanted in the channel of a silicon 
MOSFET. de Sousa has proposed a model of this noise 
as caused by dangling bonds located at the interface be- 
tween the crystalline Si and the amorphous oxide |14j . 
These defects, known as Pi, centers, are associated with 
oxygen vacancies in the oxide and have the structure 
•Si = Sis- The lone electron in the dangling bond can pro- 
vide a thermally switching paramagnetic contribution to 
the magnetic environment experienced by the donor that 
causes magnetic field noise. This noise then acts to de- 
phase the electron spin qubit defined on the phosphorus 
dopant. By modeling the Pf, center spin flips as coupled 
to tunneling two-level systems in the oxide, de Sousa has 
shown that the resulting magnetic field noise possesses 
an approximate power spectral density ^(a;) oc l/co |14j . 
Recent work by Paik et al. has provided evidence in favor 
of de Sousa's model [T5] . 

We can empirically estimate the proportionality con- 
stant, or equivalently the noise strength e (Eq. ^) by 
comparing the calculated T2 time of an electron spin 
qubit with the experimental values extracted for phos- 
phorus donors implanted in Si MOSFET devices. Donor 
electron T2 times have been shown to be several millisec- 
onds in isotopically purified silicon [121 [H] . 

We can determine an effective T2 time for a pulse se- 
quence in our model as follows. We initialize a qubit in 
the -|-1 eigenstate of ax, and apply a quantum memory 
pulse sequence repeated several times. If one measures 
only at the conclusion of each repetition of the pulse se- 
quence, the quantity (ax) will decay approximately ex- 
ponentially in time, with a time constant that we define 
to be T2. 

We have set the value of the parameter e to be that 
which yields T2 ~ 1ms under a Carr-Purcell pulse se- 
quence operated with a 1% duty cycle. As demonstrated 
in Figure [3) this condition is satisfied with e = 0.0011, 
which we have rounded to e = 10^"^. 

In the remainder of this section we explore the power of 
numerically optimized pulse sequences obtained with the 
improved multistate Markovian fluctuator, for two tar- 
get unitary operations subject to this interface-induced 
l/uj noise with an additional constant noise offset that al- 
lows for heterogeneous dephasing. The first is quantum 
memory, i.e., the preservation of coherence of an arbi- 
trary quantum state, while the second is a single-qubit 
Hadamard transformation. We find that excellent per- 
formance of GRAPE for both operations can be attained 
even in the presence of additional zero-frequency (con- 
stant) noise. 
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Figure 3: (Color online) Calculated qubit T2 times as a func- 
tion of the noise parameter e. Results are presented for qubit 
T2 with the application of a Carr-Purcell sequence with a 1% 
duty cycle. The intersecting horizontal and vertical lines in- 
dicate the value of e at which the electron T2 equals 1ms. 



Consistent with the application to experimentally ac- 
cessible phosphorus dopants silicon devices, we construct 
here pulse sequences that may be implemented by cur- 
rent signal generators and that are thus subject to lim- 
itations on the on/off ratio. Consequently, we enforce a 
50% duty cycle on the GRAPE sequences, i.e., each pulse 
is followed by a brief quiescent period and where the to- 
tal quiescent time must be at least half of the total pulse 
length, or, r/2. 



Quantum Memory 

We begin this section with a discussion of coherence 
maintaining operations, known generally as quantum 
memory gates. When designing such pulses, one must 
make pulse design decisions based on the specifics of the 
experiment in question. Consider a particular experi- 
ment which requires that coherence be maintained for a 
certain time, t. Ideally, one would design a pulse sequence 
itself having total length t, as well. In principle, such a 
pulse sequence would have much more flexibility than a 
sequence of duration t/n repeated n times. However, the 
computational effort scales exponentially in the number 
of pulses, as discussed above, and long sequences may be 
difficult to find which match or exceed the performance 
of repeated short sequences. Which choice is made will 
depend strongly on the computational resources available 
to the pulse designer. 

Here we present a numerical solution of Eq. ^ with G 
the identity matrix, total pulse sequence time t = SOr^n- 
and total number of pulses iVp = 30. The particular op- 
timal solution found under these constraints is shown in 
Fig. |4] Here r^ refers to the time required to perform a 
full TT rotation of the qubit at maximum control ampli- 
tude. This value of A^p was chosen because it was the 



largest value for which we were able to obtain results in 
a reasonable amount of computer time (see discussions 
above). To compare with these numerically optimized 
pulse sequences we construct an equivalent length Carr- 
Purcell (CP) decouphng pulse sequence, defined by 



W-~TTr — W — W — TTr 



repeated 7 times, where n^ denotes a 7r-pulse about the 
X-axis, and w denotes a quiescent period of t„, — jT^^. 
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rotations, which prevents the exact refocusing of zero- 
frequency noise, resulting in a significant rise in error as 
the zero-frequency noise magnitude \r]os\ increases. Note 
that the lack of time reversal symmetry possessed by the 
optimized pulses (unlike the Carr-Purcell sequence) re- 
sults in an asymmetry with respect to 770s, as illustrated 
in Fig. [5] 

The numerically optimized sequence shows improved 
performance relative to these Carr-Purcell decoupling se- 
quences in two respects. First, it performs better than 
the Carr-Purcell pulse sequences for zero and small |?yos|, 
due to the greater flexibility of the numerical optimiza- 
tion in developing protection against the 1/uj noise com- 
ponent. Thus, at ?7os = 0, the error obtained with the 
numerically optimized pulse sequence is 2.88 x 10^^, com- 
pared to 3.26 X 10~^ with the finite amplitude Carr- 
Purcell pulse sequence. However, a far more dramatic 
difference is the greater robustness against the magni- 
tude of zero-frequency noise. The numerically optimized 
pulse sequence is seen to show very small error over a 
broad range of 770s , attesting to the power of the numer- 
ical approach to mitigate combined decoherence effects 
deriving from very different noise sources. 
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Figure 4: (a) a;-axis control function and (b) y-axis control 
function for implementing quantum memory operations. 

Fig. [5]shows the error (defined as 1 - $/) as a function 
of the zero frequency noise rjos for a numerically gen- 
erated pulse sequence that is optimized over all values 
of rjos (red line), in addition to optimization against the 
1/cj noise. The blue and green lines show the correspond- 
ing fidelities obtained with the Carr-Purcell sequence of 
equivalent duration specified above, using finite ampli- 
tude (dot-dashed blue line) and infinite amplitude (dot- 
ted green line) pulses. 

Infinite amplitude Carr-Purcell pulses are capable of 
refocusing arbitrarily large zero-frequency noise, result- 
ing in a constant error as a function of ?7os whose value can 
be taken as a measure of the uncorrected error due to the 
I /id noise component. Unlike the ideal, infinite ampli- 
tude pulse sequence, a Carr-Purcell sequence with finite 
amplitude pulses does not allow complete Bloch sphere 
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Figure 5: (Color online) Error (1 — $/) in quantum memory 
of a qubit subject to dephasing noise with spectral density 
St oc 1/uj + 77os(5(tj) under various control pulse sequences of 
duration SOttt, shown as a function of offset noise values 770s. 
The solid red line represents the error for a sequence that 
is optimized over a range of offset noise. Green dotted line: 
error obtained with infinite-amplitude Carr-Purcell sequence. 
Blue dot-dashed line: error obtained with finite amplitude 
Carr-Purcell sequence. Also shown as the black, dashed line 
is a pulse sequence obtained without regard to robustness over 



Hadamard Gate 

Our second target operation is the Hadamard gate. 
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a common single qubit operation in quantum algorithms. 
The optimization considerations for implementing such a 
single qubit rotation with numerically optimized pulse se- 
quences are similar to those for generating sequences to 
protect the identity gate. However, in contrast to the 
situation for quantum memory, here we are interested in 
maximizing fidelity and robustness to a constant noise 
offset, rather than in the maintenance of coherence over 
a long time. Thus the optimal pulse sequences for protec- 
tion of the Hadamard gate are considerably shorter than 
the sequences derived above for protection of quantum 
memory. 

Using the same cost function as Eq. (|8| and evaluating 
^H for the H operation, we were able to obtain high- 
fidelity pulse sequences with r = 6t^ and A^p = 6. Fig. [6] 
shows the resulting pulse sequence when optimization is 
made for the case of zero offset noise, ijos = 0. 



with a single pulse sequence that is optimized only for 
rjos = shows again the enhanced robustness afforded by 
the numerical optimization approach. 
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Figure 6: Two-dimensional control function producing a high 
fidelity H rotation inT — 6n for rjos = 0. This pulse sequence 
results in a worst case fidelity of $h = 1 — 8.27 x 10~^ and 
exhibits a strong robustness to the value of constant offset 
noise rjoa. Panel (a) shows the a;-axis control function and 
panel (b) the y-axis control. 



This pulse sequence results in a worst case error for 
the Hadamard gate of 8.27 x 10~^ at rjos = 0. Thus it 
is evident that gate operations can readily be corrected 
at similar or better levels than quantum memory, using 
shorter pulse sequences. Fig. [7] shows the worst case er- 
ror for a numerically optimized pulse sequence optimized 
over a range of 770s values, as before, as a function of the 
noise offset 770^. Comparison with the results obtained 
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Figure 7: (Color online) Error (1-$h) of a _ff rotation un- 
der numerically optimized pulse sequences mitigating against 
noise with spectral density St oc l/uj + ri'^^5{Lo) as a function 
of rios- The solid, red line shows the result of optimizing with 
respect to rjos as well as the 1/uj noise. The dashed, black line 
shows the considerably less robust result of using a pulse se- 
quence that is optimized only at a single value of rj (rjos ~ 0) . 



SUMMARY 

We have expanded the techniques of (TJHIE] to develop 
a general numerical method for simulating noise sources 
deriving from a broad variety of Markovian power spec- 
tra. The method employs a new numerical approach to 
generation of the noise spectrum that can greatly reduce 
the number of noise states required to simulate a noise 
source with a given spectral density. We illustrated this 
with the example of a four state simulation of a 1/lu spec- 
trum over two decades of frequency w, which is signifi- 
cantly more efhcient than the constructive method em- 
ployed previously in Ref. 0]. This numerical representa- 
tion of Markovian noise was then used in the context of 
numerical generation of control pulse sequences to miti- 
gate the effects of this noise on a single qubit. Here we 
have extended the methods of [4J to allow control pulses 
to be performed along both x- and y-axes of the qubit, 
rather than along a single axis. Using numerical genera- 
tion of the noise allows spectral densities from different 
sources of dephasing to be combined, giving rise to con- 
siderable additional flexibility and robustness in the de- 
coherence mitigation. This was illustrated by generation 
of pulse sequences designed to minimize decoherence in 
the presence of both homogeneous dephasing character- 
ized by 1/uj spectral density and a source of heteroge- 
neous dephasing, characterized by a zero frequency noise 
offset rjos- The numerical optimization approach allows 
the pulse sequences to be simultaneously optimized with 



respect to the parameter 770s and the 1/w noise. This 
introduces an unprecedented robustness to decoherence 
mitigation with reahstic bounded amphtude controls in 
the presence of resonance frequency errors and inhomo- 
geneous broadening. In particular, the performance of 
the numerically optimized pulse sequences over a range 
of 770s values was seen to be considerably superior to the 
corresponding performance of a standard dynamical de- 
coupling pulse sequence with bounded amplitudes. 

As a demonstration of the power and flexibility of these 
numerical methods for noise mitigation, we have explic- 
itly studied the protection of quantum memory and the 
protection of the Hadamard gate. To ground the derived 
pulse sequences to a physical system, we took estimates 
of noise strength that are appropriate to the situation of 
dephasing noise acting on phosphorus donors in silicon 
and implemented the numerical optimization subject to 
realistic constraints of duty cycle and pulse amplitude 
limitations. The remarkable robustness of the optimal 
pulse sequences with respect to the constant noise offset, 
showing worst case gate errors of order 10^^ — 10^^ over 
a range of noise offsets, is encouraging for application of 
these pulse sequences to current experiments with spin 
qubits in semiconductors [16]. 
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